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The QCD phase transition at high temperature and low density 

S. Ejirif C.R. AUtonf S.J. Hands? O. KaczmarekJ' F. KarschJ' E. LaermannJ' and Ch. Schmidt*^ 
'^Department of Physics, University of Wales Swansea, Singleton Park, Swansea, SA2 8PP, U.K. 
'^Fakultat fiir Physik, Universitat Bielefeld, D-33615 Bielefeld, Germany. 

We study the thermal properties of QCD in the presence of a small quark chemical potential fi. Derivatives of 
the phase transition point with respect to fi are computed at /i = for 2 and 3 flavors of p-4 improved staggered 
fermions on a 16^ x 4 lattice. Moreover we contrast the case of isoscalar and isovector chemical potentials, quantify 
the effect of /i 7^ on the equation of state, and comment on the screening effect by dynamical quarks and the 
complex phase of the fermion determinant in QCD with /i 7^ 0. 



To understand recent heavy-ion collision exper- 
iments, theoretical study of the QCD phase tran- 
sition at high temperature and low density is im- 
portant. For instance, the interesting regime for 
RIHC is /iq/Tc ^ 0.1, where /iq = /i/a is a quark 
chemical potential. However, the Monte-Carlo 
method is not directly applicable for simulations 
at /i 7^ 0, which makes the study of finite-density 
QCD difficult; hence we usually use the reweight- 
ing method. Using the identity 



(0>(/3,M)= W(/3„.0)/W(/5o,0). 



(1) 
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,Arf(lndetM(/i)-lndetM(0))„-Sg(/3)+Sg(/3o) 



the expectation value (O) at /i 7^ is computed 
by a simulation at ^ = 0. Here M is the fermion 
matrix, Sg the gauge action, and A^f the number 
of flavors. Then, there exists a famous "sign prob- 
lem". Because detM is complex at /i 7^ 0, if the 
complex phase fluctuates rapidly, both numera- 
tor and denominator in RHS of eqn.(|l|) become 
vanishingly small. For the case of small /i, the 
complex phase can be written by the odd terms 
of the Taylor expansion of In det M [Q| . Denoting 
detM = IdetMle*", 
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* Presented by S. Ejiri. Numerical work was performed 
using the UKQCD APEmille in Swansea supported by 
PPARC grant PPA/A/S/ 1999/00026. 



The first term is N{lmtr[M-^{dM/diJ,)]ii. From 
these equations, we find explicitly that the mag- 
nitude of 9 is proportional to /x, the volume and 
N{. Moreover the first term of 9 can be com- 
puted by the noise method and has a tendency 
that the phase fluctuation becomes larger as the 
quark mass or T decreases. Roughly speaking, 
the sign problem happens when the fluctuation 
of 9 becomes larger than 0(7r/2). Thus the sign 
problem is not serious for the range of small /it, 
but that region becomes narrower and narrower 
as the volume increases, which suggests that in 
the thermodynamic limit, since the region where 
the sign problem is manageable decreases to zero 
size, the only way of accessing the /i 7^ region is 
via a Taylor expansion since this involves calcu- 
lating quantities (i.e. derivatives) only at /i = 0. 

Last year, we proposed a general formulation to 
compute the derivatives of physical quantities . 
We perform a Taylor expansion for In det M and 
fermionic operators in eqn.(0) and neglect higher 
order terms of fi. Then the resulting expectation 
value contains an error of higher order in /i but 
which docs not affect the calculation of deriva- 
tives of lower order than the neglected terms. In 
this report, we summarize the results obtained by 
this method. We perform simulations on a 16^ x 4 
lattice using a combination of the Symanzik im- 
proved gauge and 2 flavors of the p4-improved 
staggered fermion actions Details are writ- 
ten in Ref. m. We also comment on the iVf de- 
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Figure 1. Polyakovloop susceptibility atm — 0.2. 



pendence, in the final section, by performing an 
additional 3 flavor simulation. 

Phase transition line 

First of all, we discuss the phase transition line 
in the (T, fi) plane. Because the first derivative 
is expected to be zero from the symmetry un- 
der exchange of fi and —fi, we calculate second 
derivative of Tc with respect to /z. In Fig. |l|, we 
plot Polyakov loop susceptibility as a function of 
/3 for ^ = 0, ±0.05 and ±0.1 at m = 0.2. This 
calculation contains errors at 0{fj,^). From this 
figure, we find that the peak position of the sus- 
ceptibility moves left as ^ increases, which means 
that f3c or Tc decreases as fj, increases. Assuming 
the peak position is at Pc, we determine the sec- 
ond derivative of f3c- Combining with results for 
the chiral susceptibility at m = 0.1 and 0.2, we 
obtain d'^Pc/d^l'^ « -1.1 with 30 - 50% error and 
any quark mass dependence of d^/3c/d/x^ is not 
visible within the accuracy of our calculation. 

The second derivative of Tc is given by 



d^Tc 
d^l 



1 d2/3e 



N?Tc 



d/j.^ 



d(3 
da 



(3) 



with a(d/3/da) obtained from the string tension 
data in (|. We then find Tdd^Tc/d^l) « -0.14. 
We sketch the phase transition line from the cur- 
vature with 50% error in Fig At the relevant 
point for RHIC, this shift of Tc is very small from 
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Figure 2. Sketch of the phase diagram. The di- 
amond symbol is the critical point obtained by 
Fodor and Katz [Q. Dotted line is upper bound 
of the fit range to determine the curvature. 



that at ^ = and the result is roughly consistent 
with those obtained by other groups 

Screening effect by dynamical quarks 

Under the presence of chemical potential, the 
time reversal symmetry is broken. By the asym- 
metry, an interesting property is shown in the 
measurement of the Polyakov loop, an external 
quark current running in the positive time direc- 
tion. Negative chemical potential induces the dy- 
namical generation of anti-quarks, which in con- 
trast to quarks can completely screen an exter- 
nal color triplet current. Thus the free energy of 
a single quark is reduced, especially in the con- 
finement phase, and the singularity at the phase 
transition point is weakened, since the long range 
fiuctuation is screened. This effect can be seen in 
Fig. |l], where we denote the Polyakov loop sus- 
ceptibility Xi at /X < by dot-dot-dash and dot- 
dash-dash lines. We see that the peak height of 
XL becomes smaller for /i < corresponding to 
a weaker singularity, while the peak position is 
almost the same between positive and negative 
/i. Simultaneously, we observe that the Polyakov 
loop at /X < is larger than that at ^ > 0, suggest- 
ing the free energy of a single quark is reduced. 
Note that this property cannot be seen in 2-color 
QCD where systems at fj. and — /i are identical. 
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Equation of state 

Next we discuss ttie /u-dependence of the equa- 
tion of state which describes the energy density e 
and pressure p. If we employ the integral method 
based on the homogeneity of the system, we ob- 
tain p = where / = -(T/T/)lnZ. Then 
derivatives of p with respect to fi are related to the 
quark number density Uq and the singlet quark 
number susceptibility xs = duq/dfiq 



d{p/T^) 

dfiq 



y4 ' 



dHp/T*) 



Xs 
2^4 ■ 



(4) 



The quark number density is zero at ^ = so the 
leading correction is 0{fi^). Moreover the second 
derivative of e can also be estimated by 



a'(e-3p)/T^ 



1 dxs f I da 



dp \adp 



(5) 



Here we neglect a{dm/da), an approximation 
which is valid in the chiral limit. Wc obtain 
T^^^{p/T^)/^^il « 0.69 and T^d^{e/T'^)/dtil « 
10.6 at Pc for TO = 0.1. The discrepancy oi p/T'^ 
(e/T^) in the RHIC regime fiq/T^ - 0.1 from its 
value at = is about 0.0035 (0.05). This is a 
1% effect, and hence quite small. 

From the second derivatives of p and e with re- 
spect to fi together with the derivatives with re- 
spect to T, we calculate the line of constant pres- 
sure and energy density. We find that the slope 
of the constant pressure (energy density) line is 
T{dT/d{nl)) ~ -0.11(-0.09). As the slope of 
in fi^ is Tc(dTc/d(/i2)) « -0.07, this result sug- 
gests that the line of constant pressure or energy 
density is parallel with the phase transition line. 

Iso-vector chemical potential 

If instead we were to impose an isovector chem- 
ical potential fij having opposite sign for u and d 
quarks , then the quark determinant would be- 
come real and positive, enabling simulations using 
standard Monte-Carlo methods In the frame- 
work of the Taylor expansion, terms even in fi 
are identical for both u and d quarks, but odd 
terms cancel for the case /x/ 7^ 0. We analyzed the 
transition point PcifJ-i) and do not observe signif- 
icant difference for Pc between fi and fj,i in the 
region of small /i, which is different from a naive 



expectation that the phase transition line for /x/ 
runs toward the T = onset threshold of a pion 
condensate at a critical fijo — rnps/^ < ttln/S, 
and hence the curvature for ^/ larger than for 
isoscalar /x. However the quark mass we used is 
still large and TOtt is not so small. Simulations at 
small TO are necessary to check the naive picture. 

A^f dependence 

Finally, we comment on the difference between 
N[ — 2 and 3. We performed an additional sim- 
ulation of A^f = 3. The preliminary results of the 
curvature at to = 0.1 are Tdd'^Tc/dfiq) « -0.15 
and —0.11 for the cases fiud — ^J's — ^J'q and 
Mud — /^q,Ms — respectively. There is no sig- 
nificant difference from Nf = 2. Therefore we 
expect that our result for A't = 2 is not so dif- 
ferent from real QCD lying between Ni — 2 and 
3. 

The most interesting point for A^f = 3 is the ex- 
istence of a critical quark mass nic on the 11 — 
axis between a first order phase transition at 
small TO and a crossover at large m. We also ex- 
pect such a critical point at /i ^ even in the re- 
gion of large m. Hence to investigate the relation 
between these critical points is quite important. 
The first preliminary result is reported by 
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